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ABSTRACT 

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) - 
which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact 
objects in supernova remnants (CCOs) - represent a key element in understanding the neutron star phenomenology. We report the 
results of an observational campaign to study the properties of the source 2XMM J104608.7-594306, a newly discovered thermally 
emitting INS. The evolutionary state of the neutron star is investigated by means of deep dedicated observations obtained with the 
XMM-Newton Observatory, the ESO Very Large Telescope, as well as publicly available y-ray data from the Fermi Space Telescope 
and the AGILE Mission. The observations confirm previous expectations and reveal a unique type of object. The source, which is 
likely within the Carina Nebula (JVh = 2.6 x 10 21 cm -2 ), has a spectrum that is both thermal and soft, with kT m = 135 eV. Non-thermal 
(magnetospheric) emission is not detected down to 1% (3cr, 0.1 - 12keV) of the source luminosity. Significant deviations (absorption 
features) from a simple blackbody model are identified in the spectrum of the source around energies 0.6 keV and 1.35 keV. While 
the former deviation is likely related to a local oxygen overabundance in the Carina Nebula, the latter can only be accounted for 
by an additional spectral component, which is modelled as a Gaussian line in absorption with EW = 91 eV and a = 0. 14keV 
(lcr). Furthermore, the optical counterpart is fainter than ray = 27 (2<x) and no y-ray emission is significantly detected by either 
the Fermi or AGILE missions. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only 
RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, 
P = 18.6 ms, at the 4cr confidence level. We interpret these new results in the light of the observed properties of the currently known 
neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of 
these scenarios can satisfactorily explain the collective properties of 2XMM J104608. 7-594306, although it may be related to the still 
poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AOl 1), will 
help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin 
down. 

Key words, stars: neutron - X-rays: individuals: 2XMM J104608.7-594306, 1RXS J141256.0+792204, RRAT J1819-1458 



1. Introduction 

A major outcome of the ROSAT mission was the discovery of 
a group of seven radio-quiet thermally emitting isolated neutron 
stars (INSs), which were originally identified serendipitously as 
soft, bright X-ray sources with no obvious optical counterparts. 
They share a rather well-defined set of properties that have never 
been encountered together in previously known classes of INSs, 
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and have been nicknamed the "Magnificent Seven" (M7, for 
short; see Haberl 2007; Kaplan 2008; Turolla 2009, for reviews). 

Relative to standard (rotation-powered) radio pulsars, these 
INSs rotate slower (P ~ 10 s), have thermal X-ray luminosi- 
ties higher than their spin-down power, and stronger magnetic 
fields (B > 10 13 G, as inferred from timing measurements as 
well as from the broad absorption lines in their X-ray spectra; see 
Kaplan & van Kerkwijk 201 1, and references therein). The M7 
display neither persistent nor transient radio emission to a rather 
sensitive limiting flux (Kondratiev et al. 2009) and are unasso- 
ciated with supernova remnants. Proper motion studies (Motch 
et al. 2009, and references therein) show that they are cooling, 
middle-aged neutron stars (ages 10 5 to 10 6 yr), probably born 
in the nearby OB associations of the Gould Belt. Their prox- 
imity (few hundred parsecs; Posselt et al. 2007; van Kerkwijk 
& Kaplan 2007) and the combination of strong thermal radi- 
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ation and absence of significant magnetospheric activity make 
them ideal targets for testing surface emission models, deriving 
radii and constraining the equation of state of neutron star mat- 
ter. Unfortunately, the current lack of understanding of the sur- 
face composition, magnetic field, and temperature distributions 
have limited any definite conclusion so far (see e.g. Kaplan et al. 
2011). 

Growing evidence relates the M7 to other peculiar groups 
of INSs, in particular the magnetar candidates, anomalous X- 
ray pulsars, and soft y-ray repeaters (AXPs and SGRs, see e.g. 
Mereghetti 2008, for a review), and the rotating radio transients 
(RRATs, McLaughlin et al. 2006; Keane & McLaughlin 2011). 
The long spin periods, bright thermal emission, and intense mag- 
netic fields suggest that some of the M7 might have evolved from 
the younger and more energetic magnetar objects (e.g. Heyl & 
Kulkarni 1998; Popov et al. 2010). Crustal heating by means of 
magnetic field decay seems to play an important role in the ther- 
mal evolution of neutron stars with B>5xl0 13 G, since a cor- 
relation between B and kT is observed (e.g. Aguilera et al. 2008; 
Pons et al. 2009; Kaplan & van Kerkwijk 2009). On the other 
hand, the discovery of RRATs is intriguing since these sources 
have so far manifested themselves in a variety of ways. While 
many, in spite of their transient nature, display timing proper- 
ties that cannot be distinguished from those of the bulk of nor- 
mal radio pulsars, several RRATs with longer spin periods and 
higher magnetic fields occupy a region of the P — P diagram that 
is devoid of normal pulsars and close to that populated by the 
M7. Moreover, the most active source among the known sam- 
ple is the highly magnetized RRAT J1819-1458, the only one 
that has so far been detected in X-rays (Reynolds et al. 2006; 
McLaughlin et al. 2007; Kaplan et al. 2009) 1 . Very interestingly, 
the X-ray source was found to exhibit a spectrum that although 
hotter and more absorbed is remarkably similar to those of the 
M7. Unusual timing behaviour following glitches detected in the 
radio indicates that this source might also have evolved from a 
magnetar (Lyne et al. 2009). 

It is remarkable that a group of very similar sources, display- 
ing at the same time unique properties that are extremely differ- 
ent from those of standard radio pulsars, are all detected in the 
very local solar vicinity. Taking into account that the M7 repre- 
sent about half of all young (younger than 3 Myr) INSs known 
within ~ 1 kpc (Popov et al. 2003), they could be just the tip of 
an iceberg of a largely hidden population of stellar remnants. 
Both the poorly constrained populations of RRATs and ther- 
mally emitting INSs are indeed estimated at present to outnum- 
ber ordinary radio pulsars by factors of one to three (e.g. Keane 
et al. 2010), whereas the latter alone are sufficient to account 
for the total number of past Galactic core-collapse supernovae. 
Therefore, to avoid "superpopulating" the Galaxy with unrelated 
neutron stars, it is necessary to invoke links between the several 
subgroups, which may be either selection effects (i.e. the same 
object appears differently to the observer due to viewing biases) 
or actual evolutionary relations. 

1.1. The source 2XMM J1 04608. 7-594306 

Radio-quiet INSs are extremely elusive and difficult to de- 
tect. For several years, considerable efforts have been devoted 
to discovering new sources similar to the M7 (e.g. Rutledge 
et al. 2003; Chieregato et al. 2005; Agueros et al. 2006; 



1 Of the more than 60 RRATs known to date, ~ 20 have precise posi- 
tion determinations and six have been investigated in X-rays; see Keane 
et al. (201 1) and Keane & McLaughlin (201 1) for reviews. 



Turner et al. 2010; Agueros et al. 2011). These involve cross- 
correlating X-ray sources - often with large positional errors 
- with a large number of catalogued objects at other wave- 
lengths, and dealing appropriately with the level of sample con- 
tamination by other classes of X-ray emitters. In spite of these 
searches, only two radio-quiet thermally emitting INSs have 
been positively identified to date since the ROSAT era, namely 
Calvera (1RXS J141256.0+792204) and 2XMM J104608.7- 
594306 (hereafter J1046). 

Calvera is a relatively bright ROSAT source, which was orig- 
inally identified as a likely compact object based on its large X- 
ray-to-optical flux ratio (Rutledge et al. 2008). The distance to 
the source is poorly constrained. Its interpretation as a cooling 
INS (with a similar size of emission radius as the M7) has always 
met with difficulties, since its high Galactic latitude and hot tem- 
perature (kT ~ 200 eV) would either require a very high spatial 
velocity to explain its current position well above the Galactic 
plane or point to a non-standard cooling or re-heating process. 
Alternatively, if Calvera's progenitor were a high-velocity run- 
away star, it might have been born away from the Galactic plane 
(Posselt et al. 2008). Recent multiwavelength investigations by 
Zane et al. (201 1) confirmed the radio-quiet nature of the source 
and most notably revealed the neutron star's very short spin pe- 
riod of P ~ 59 ms, which was later also apparently detected in 
y-rays with data from the Fermi Large Area Telescope (LAT, 
Atwood et al. 2009). The nature of Calvera however remains un- 
clear, especially after Halpern (201 1) called the y-ray detection 
into question, leaving the spin-down rate of the source far more 
poorly constrained (see Sect. 5). 

We previously reported the results of a programme to iden- 
tify new thermally emitting INSs in the XMM-Newton cata- 
logue of serendipitous X-ray sources (Pires et al. 2009a). The 
X-ray brightest INS candidates resulting from our search have 
been the scope of succesful proposals carried out over the past 
four years in the optical with the 8.1m European Southern 
Observatory Very Large Telescope (ESO-VLT), the two 8.4m 
Large Binocular Telescope (LBT), and the 4.1m Southern 
Observatory for Astrophysical Research (SOAR). These deep 
optical investigations permitted the discovery of a new thermally 
emitting INS, J1046 (Pires et al. 2009b, hereafter Paper I) . 

As expected for a cooling INS, J1046 is characterised by a 
soft blackbody-like emission, stable observed flux over a long 
timescale, a very high X-ray-to-optical flux ratio, and no radio 
counterpart in the Parkes Multibeam Pulsar Survey (which has a 
sensitivity of 0.14mJy for a canonical pulsar). The column den- 
sity is consistent with the source being located within the Carina 
Nebula, at a distance of 2 kpc (Paper I). This giant HII region 
is indeed likely to harbour unidentified cooling neutron stars, 
which are the faint analogues of those produced in other close- 
by regions of intensive star formation (Posselt et al. 2008, see 
their Fig. 7). 

To investigate its properties and evolutionary state, we have 
ongoing dedicated projects in the optical and X-rays to observe 
this source. We report here the results of an observational cam- 
paign with XMM-Newton and ESO-VLT. The paper is struc- 
tured as follows: in Sects. 2 and 3, we describe the X-ray and 
optical follow-up observations, together with detailed analysis 
and results. A search for pulsations in Fermi-LAT and AGILE 
(Tavani et al. 2009) data is described in Sect. 4. The summary 
and discussion of the results are in Sect. 5. Finally, our conclu- 
sions are given in Sect. 6. 
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Table 1. Instrumental configuration and duration of the EPIC scien- 
tific exposures of the XMM-Newton observation of 2XMM J104608.7- 
594306 



Instr. 


Start Time 


Mode 


Durat. 




(UTC) 




(s) 


pn 


2010- 12-06T00: 14: 10.0 


Prime small window 


90,471 


MOS1 


2010-12-06T10:01:39.0 


Prime partial W2 


55,120 


MOS2 


2010-12-06T10:01:45.0 


Prime partial W2 


55,129 



Notes. The cameras were operated in imaging mode and the thin filter 
was used (obsid 0650840101). 



2. XMM-Newton observations 

Being located at an angular distance of ~ 8.5' from the well- 
studied binary system Eta Carinae (77 Car), J 1046 has been 
detected on many occasions by XMM-Newton and Chandra. 
Unfortunately, in most observations, conditions were far from 
optimal, given the source's large off-axis angles, ~ 9' on aver- 
age, and short effective exposure times, usually less than 15 ks 
(see Paper I, for details) 2 . We were granted XMM-Newton time 
(90 ks, A09) to help us to characterize the source spectral energy 
distribution and search for pulsations. 

2.1. Observations and data reduction 

Observations were carried out on 2010 December 6, for a to- 
tal exposure time of 90.917ks. All instruments on-board XMM- 
Newton (Jansen et al. 2001) were active during the observation. 
Table 1 contains information on the scientific exposures and in- 
strumental setup of the EPIC pn (Struder et al. 2001) and MOS 
(Turner et al. 2001) detectors. 

The EPIC cameras operated in small window (SW) mode 
with thin filters. In this configuration, pn and MOS provide a 
time resolution of 6 ms and 0.3 s, respectively; owing to the more 
frequent readout, the nominal exposures are shorter by factors of 
0.71 and 0.975 (corresponding to deadtimes of 29% and 2.5%). 
For pn, only the 63 x 64 pixels (4.3' x 4.3') of CCD 4 are read 
out. For MOS, the central 100 x 100 pixel area (1.8' x 1.8') of 
the inner CCD operates with a higher time resolution, while the 
outer chips are read out normally. 

The MOS cameras experienced sporadic full scientific 
buffers during the first 35 ks of the observation, alternating very 
short scientific exposures, of typical durations of shorter than 
half a kilosecond, with periods where the cameras operated in 
the so-called "counting mode", lasting between 2.6 ks and 6.8 ks. 
In this mode, no transmission of information for individual X- 
ray events occurs. As a result, the effective observing time is 
considerably shorter for MOS. The reason for the abnormal op- 
eration of the MOS detectors is the presence of the X-ray bright 
77 Car in one of the peripheral MOS CCDs, which are used to 
monitor the radiation level during the observation under the as- 
sumption that no bright sources are present. The MOS expo- 
sures were unfortunately interrupted and restarted several times 
until the origin of the high count-rate level could be identified 
(R. Gonzalez-Riestra, private communication). 

Data reduction was performed with SAS 11 (xmmsas- 
20110223-1801) by applying standard procedure and using the 
latest calibration files. MOS exposures U008 and pn S003 were 

2 In spite of the ten-year-long timebase, a proper motion study in X- 
rays is impractical due to the source's estimated distance (d ~ 2kpc) 
and to the large off-axis angles and corresponding positional errors in 
the observations before A09. 



Table 2. Source coordinates with lcr errors 



Camera 


KA 


DcC 


<x 




(h m s) 


(d m s) 




pn 


10 46 08.500 


-59 43 05.179 


0.2 


MOS1 


10 46 08.536 


-59 43 04.863 


0.3 


MOS2 


10 46 08.485 


-59 43 04.451 


0.3 


EPIC 


10 46 08.503 


-59 43 05.143 


0.2 


EPIC T 


10 46 08.756 


-59 43 05.523 


0.4 


Chandra 


10 46 08.719 


-59 43 06.480 


< 1 


XMM-Newton 


10 46 08.730 


-59 43 06.300 


0.4 



Notes. We also list the position of the source as previously determined 
in archival Chandra and XMM-Newton data. f Corrected position after 
cross-correlating the EPIC source list with optical objects in the field 
(see text). 



processed using the EPIC meta tasks emchain and epchain, re- 
spectively, applying default corrections. We ensured that the pn 
event files were clean of unrecognised time jumps (i.e. those un- 
corrected by standard SAS processing). 

No background flares were registered during the 90 ks (last 
55 ks) of the pn (MOS) observations. The effective observing 
times were therefore only reduced by the cameras' livetime in 
SW mode: 63.25 ks for pn and 53.4ks for MOS. We note that 
although optical loading (due to unblocked optical photons from 
77 Car and other bright OB stars in the field-of-view) is obvi- 
ously present at r\ Car's position in the MOS images, out-of-time 
events or charge transfer efficiency (CTE) alteration are unim- 
portant at the position of the target at aimpoint, since 77 Car is 
not located on the same CCD as J 1046 or along the readout di- 
rection (M. Ehle, private communication). 

For the analysis, we filtered the event lists to retain the pre- 
defined photon patterns with the highest quality energy calibra- 
tion - single and double events for pn (pattern < 4) and single, 
double, triple, and quadruple for MOS (pattern < 12) - as well 
as to exclude bad CCD pixels and columns. Source photons were 
extracted from circular regions of radius 18" (unless otherwise 
noted) centred on the position of the source (as given by emlde- 
tect for each camera; see Sect. 2.2). Background circular regions 
of size 50" (pn) and 25" (MOS) were defined off-source, on 
the same CCD as the target. Owing to the frequent readout of 
the pn in SW mode, detector noise dominates the low-energy 
count distribution. We therefore restricted the analysis of pn data 
to photons with energies between 0.3 keV and 2 keV (see how- 
ever Sect. 2.3). For MOS, the energy band for the analysis is 
0.15 -2keV. 

The detected source count rates are 7.74(14), 1.66(6), and 
1.90(6) (xl0~ 2 s" 1 ; pn, MOS1, and MOS2); these closely agree 
with the expectations of simulations (considering the source pa- 
rameters as in Paper I). Using the SAS task epiclccorr, we created 
lightcurves for J 1046, which were corrected forbad pixels, dead- 
time, exposure, and background counts. The statistics for these 
corrected lightcurves, binned into 900 s and 550 s (pn and MOS) 
intervals, show the 3cr upper limits for the rms fractional varia- 
tion of 0.12, 0.30, and 0.25 for pn, MOS1, and MOS2, respec- 
tively. On the basis of the same lightcurves, the reduced xl as- 
suming a constant flux is 1.007 (pn, 100 d.o.f.), 0.979, and 1.067 
(MOS1 and MOS2; 99 d.o.f.), corresponding to null-hypothesis 
probabilities of 46%, 5 1 %, and 28%. 
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Table 3. Parameters of 2XMM J104608.7-594306, as extracted from the XMM-Newton AQ9 observations 



Parameter 


pn 


MOS1 


MOS2 


EPIC 


Counts 


5594 + 91 


1024 ±41 


1133 ±40 


7912 ± 108 


Counts (0.2 - 0.5 keV) 


1236 + 41 


154 ± 15 


161 + 15 


1258 + 42 


Counts (0.5-1.0keV) 


3437 + 71 


611 ± 32 


627 ± 29 


3555 ±71 


Counts (1.0- 2.0 keV) 


910 ±37 


257 ± 21 


343 ± 22 


921 ± 37 


Counts (2.0 - 4.5 keV) 


0±7 


2±4 


1.4 ±2.5 


2±5 


Counts (4.5- 12 keV) 


11 ± 12 


0±4 


0±5 


14 ± 12 


Detection likelihood 


6522 


968 


1519 


9436 


Flux (xl0~ 13 ergs -1 cirT 2 ) 


1.17(4) 


1.09(7) 


1.17(7) 


1.174(29) 


Rate(xl0- 2 s -1 ) 


9.03(15) 


2.00(8) 


2.17(8) 


13.43(18) 


I (degrees) 


287.73 


287.73 


287.73 


287.73 


b (degrees) 


-0.59591 


-0.59579 


-0.59574 


-0.59587 


HRi 


+0.471 ±0.015 


+0.60 + 0.04 


+0.59 + 0.03 


+0.477 ± 0.015 


HR 2 


-0.581 ±0.015 


-0.41 ±0.04 


-0.29 + 0.04 


-0.588 ±0.015 


HR 3 


-1.000 ±0.016 


-0.985 + 0.028 


-0.992 ±0.014 


-0.997 ±0.011 



Notes. Counts, fluxes, and rates are given in the total XMM-Newton energy band (0.2 - 12keV). 



2.2. Source detection 

To measure the position of J 1046, we created images in the five 
XMM-Newton pre-defined energy bands, for each EPIC cam- 
era, and analysed them (individually and simultaneously) with 
the source-detection meta task edetecLchain. In Table 2, we 
list the source equatorial coordinates in each camera together 
with previous determinations based on Chandra and XMM- 
Newton (obsid 031 1990101, with the highest signal-to-noise ra- 
tio) archival data. 

The position of the source was found to be offset rela- 
tive to previous determinations. By cross-correlating the list of 
(pipeline-processed) EPIC X-ray sources with those of optical 
(GSC2.3) objects lying within 17' from J1046, we found offsets 
of -1.8" in RA and 0.9" in DEC, as well as a slight rotational 
offset of 0.013°, using the SAS task eposcorr and a number of 59 
matches. The corrected source position, updated accordingly, is 
consistent with previous determinations (Table 2). We calculated 
the lcr errors in the updated position by quadratically adding a 
systematic uncertainty of 0.35" to the nominal errors computed 
by emldetect, which is the procedure adopted by the catalogue 
pipeline of XMM-Newton, based on correlations with optical 
catalogues. 

Source parameters based on a maximum likelihood fitting 
are given in Table 3. The hardness ratios and corresponding er- 
rors were computed according to 



where C, are counts in a given energy band and cr c are the cor- 
responding errors. Owing to the low counts detected in energy 
bands 4 and 5 (above 2keV), the errors in HR4 are too large; 
hence we only list in Table 3 the values of hardness ratios that 
can be constrained by the errors (HRi^)- 

2.3. Timing analysis 

Previous X-ray observations of J 1046 did not reveal pulsations 
to a poorly constraining 30% upper limit (3cr confidence level), 
in the 0.15 - 100 s period range (Paper I). Simulations showed 
that for an exposure time of 90 ks the EPIC cameras operating 
with the thin filter can collect enough photons to significantly 



(> 4cr) detect pulsations from J1046 of periods in the range 
P = 0.6 - 1000 s, for pulsed fractions higher than ~ 10 - 15%. 
The time resolution of the EPIC cameras in SW mode allows 
us to search for pulsations with frequencies as high as 1 .7 Hz - 
similar to those observed in both RRATs and the M7. The higher 
resolution of pn in SW mode also allows us to search for higher 
frequencies, up to 88 Hz, albeit with lower sensitivity. 

For the timing analysis, we considered events with pattern 
12 or lower. The times-of-arrival of the pn/MOS photons were 
converted from the local satellite to the solar system barycen- 
tric frame using the SAS task barycen and the source coordinates 
(Table 3). A Z\ (Rayleigh) test (Buccheri et al. 1983) was used, 
which is appropriate when searching for smooth pulsations sim- 
ilar to those in the X-ray emission of thermally emitting INSs (in 
particular, four of the M7 show pulsed fractions, pt, lower than 
10%). 

We analysed the EPIC cameras together in the P = 0.6 - 
10000 s period range, by testing different energy bands and ex- 
traction region radii. Given the high noise of the pn camera 
at low energies and the energy-dependent signal-to-noise ratio 
(S IN), we performed searches testing energy intervals where we 
aimed to achieve the optimal compromise between sensitivity 
and source signal over background counts. More specifically, 
we defined energy bands where the S/N was both the highest, 
S/N > 5 (0.4 - 1.3 keV, ~ 4260 - 5890 pn events), and mod- 
erate, S/N > 3 (0.33 - 1.5 keV, ~ 4640 - 6500 pn events); we 
also extensively tested other energy configurations by varying 
the lower (e?) and upper (e„) ends of the energy band, e; be- 
tween 0.2 keV and 1 keV and e u between 1.9keV and 2.5 keV. 
The radii of the extraction regions were varied between 10" and 
25" around the position of the source, as determined for each 
camera with emldetect (Table 2). To search for higher frequency 
pulsations, with periods as short as P — 0.0114 s, we restricted 
the analysis to the pn data, and tested different parameters (en- 
ergy bands, extraction radii, and photon patterns). The adopted 
step in frequency was 1 /iHz when analysing EPIC and twice 
this value when running the test on pn, typically for a frequency 
range ~ 50 times broader. These values oversample the expected 
width of the peaks by a factor of nearly ten, warranting that a 
peak corresponding to a periodic signal is not missed. The num- 
ber of statistically independent trials depends on the chosen fre- 
quency range of the search and the exposure time of the observa- 
tion; typically it is ~ 10 5 for EPIC searches and between ~ 10 5 
and 6 x 10 6 when screening pn data. 



4 



A. M. Pires et al.: The peculiar isolated neutron star in Carina 




0.01 



- 1 - 

0.1 



10 



— I — 

100 




1000 



0.01 



0.1 



- 1 - 

10 



100 



1000 



P(s) 



P(s) 



Fig. 1. Results of theZf search (pn data, P = 0.011 - 1000 s). The frequency range is Av = 87.72 Hz, the energy band is 0.36- 2.2 keV, and the size 
of the extraction region is 18.85" (5260 counts). Dashed horizontal lines show confidence levels of from lcr to 4<x for the detection of a periodic 
signal, given the frequency range and length of the observation. A periodic signal at P, ~ 18.6 ms is detected at 4<x (left). The periodogram on the 
right shows the same analysis conducted for a background region (size 52" and 5322 total counts; see text). 



In the frequency range where we are most sensitive to pulsa- 
tions, namely periods P = 0.6 - 1 0000 s (when the EPIC cameras 
are analysed together), no significant periodicities were found. 
The 3cr upper limit to the source pulsed fraction in this period 
range is 9.6%, considering the 0.3 - 2.5 keV energy band. 

Very interestingly, extending the search to higher frequencies 
(pn data only), we found evidence of a periodicity at P = P„ = 
18.640787(6)ms. The peak at frequency v„ = 53.645804(16) Hz 
has a maximum power of ~ 50 and corresponds to a detec- 
tion at Act (see below). The resulting periodogram can be seen in 
Fig. 1. 

Across the entire frequency range and parameter space cov- 
ered by our timing analysis, no other frequency had a power as 
high as Zj;. The peak was also the strongest for searches con- 
ducted with different photon patterns and frequency resolution. 
We found that for low frequency resolution searches, with steps 
of 8 - 10 iiliz, the exact frequency was missed and the power 
(although the highest in the search) is lower, ~ 33. Best re- 
sults were found when only single and double pn events were 
included in the search. 

To investigate whether this signal is associated with un- 
known instrumental effects, we performed the same analysis on 
photons extracted from a few background regions, with roughly 
the same total number of counts as collected for J 1046. We found 
that the maximum peaks (which occur at random frequencies) 
are always below the lcr confidence level, considering the to- 
tal frequency range of the analysis (P = 0.0114 - 1000 s and 
Av = 87.72 Hz; see again Fig 1). The corresponding power of 7? x 
statistics at P = P„ is also very low in the background regions, 
^*bkg ~ Moreover, we looked for significant peaks around 
P = P» in most recent observations of other INSs performed in 
SW mode - the 2009 observation of Calvera (obsid 0601 180201) 
and the 2006 one of RRAT J1819-1458 (obsid 0406450201). 



Peaks detected at our reported frequency are weak in both cases, 



*, Calvera 



0.7andZ t 2 RRAT 



5, as shown in Fig. 2. 



Although the highest discovered in our search, the power 
(significance) of the Z 2 { statistics at P* is sensitive both to the 
choice of energy band and the size of the extraction region 
(Fig. 3). In general, we found that potential pulsations from 
the source are likely to be smeared out by noise when the en- 
ergy band is not restricted to below ~ 0.3 keV, where the back- 
ground of the pn camera dominates. In particular, for searches 
with e; < 0.25 keV, the peak at P = P, is below the 2cr con- 
fidence level (Fig. 3; left). On the other hand, when low energy 
photons are discarded and the lower end of the energy band is not 
too restrictive (i.e. e; is between 0.3 keV and 0.6 keV), the signal 
is always detected above the 3<x confidence level. The loss of 
sensitivity detected when e; is above 0.6 keV is due to the re- 
duced number of counts (between 950 and 3750 counts in the 
1 - 2.2 keV and 0.6 - 2.2 keV energy bands, respectively, com- 
pared to the more than 5000 counts in broader energy bands). 
The signal is not as sensitive to the upper end of the energy 
band, e„, and is strongest at 2.2 keV (Z 2 =49-50 for e„ between 
1.9 keV and 2.5 keV; e t = 0.36 keV). We note that by restricting 
ourselves to high-S /N energy bands (0.4- 1.3keV, S /N > 5 and 
0.33 - 1.5 keV, S /N > 3), we obtain results that are no more 
reliable than more sensitive searches (i.e. with broader energy 
bands). 

Our timing analysis carried out using photons extracted from 
circular regions with different radii revealed that the Z 2 power 
varies steeply between 17" and 20", reaching its maximum at 
18.85" (Fig. 3; right). For the 0.36 - 2.2keV energy band, the 
signal is always above 3cr for radii between 16" and 20". Other 
energy bands also show an optimal extraction radius at the same 
value, 18.85". 



5 



A. M. Pires et al.: The peculiar isolated neutron star in Carina 



Calvera (8757 counts) 

Z 2 (Pj1046) = 0.7 



1 — 

-10 



10 



P-PJ1046 (10" 7 S) 



RRAT J1819 (2968 counts) 

Z 2 (Pj1046) = 5 



1 — 

-10 



P-PJ1046 (10" 7 S) 





P-PJ1046 (10" 7 S) 



p-Pjio46 (icr 7 S) 



Fig.2. Our Z\ analysis around the periodicity at P, ~ 18.6ms. The frequency range is 5.3 x 1CT 2 Hz and the step in frequency is 0.1//Hz; the 
energy band is 0.36 - 2.2 keV and the size of the extraction region is 18.85". We present the results of the same analysis for different data sets 
observed in SW mode: at top left, we show the periodogram for Calvera (obsid 0601180201; 8757 total counts); at top right, we show that for 
RRAT J1819-1458 (obsid 0406450201; 2968 total counts). Below are the periodograms for our data: at the left we provide results for events 
extracted from a background region (5322 total counts for an extraction region of size 52") and on the right, for J1046 (total of 5260 counts). 



For the best choices of energy band and extraction region 
radius as discussed above, the significance of the pulsed signal 
is 4a-, for Av = 87.72 Hz and N = 5.55 x 10 6 independent trials 
(the probability of the measured periodicity being spurious is 
~ 6.334 x 10~ 5 ). The corresponding pulsed fraction, defined as 



.1/2 



Pi = (iZ^^/Nph) x 100%, where N p h are the total number of 
counts in the search, is ~ 14%. A pn lightcurve, folded at P, and 
corrected for background counts and other effects, can be seen 



in Fig. 4. The pulse profile, as expected from the Z 2 search, is 
fairly sinusoidal. 



2.4. Spectral analysis 

Our analysis of archival data that serendipitously detected J 1046 
(Paper I) revealed a soft spectrum characterised by blackbody- 
like emission of an average temperature kT = 117 + 14 eV, hy- 
drogen column density 7Y H = (3.5 + 1.1) x 10 21 crrT 2 , and stable 
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Fig. 3. Left. Significance of the Z\ statistics at P = P, ~ 18.6 ms, as a function of the lower end of the energy band, e ; - 2.2 keV. The size of the 
extraction region is 18.85" . Right. Significance of the Z\ statistics as a function of the size of the extraction region, for energy band 0.36 - 2.2 keV. 
The confidence levels are given relative to the entire frequency range of the timing analysis. 
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Fig. 4. Background- subtracted light curve for pn, folded at the spin pe- 
riod P, ~ 18.6 ms. The energy band is 0.35 - 2 keV and the size of the 
extraction region is 18.85". 



observed flux of f x = 1.03(6) x 10 13 ergs 1 cm 2 (errors are at 
3o-). 

The analysis of A09 data is based on source spectra ex- 
tracted from a region of size 19", together with respective re- 



sponse matrices and ancillary files for each of the EPIC cameras 
(created using the SAS tools rmfgen and arfgen). Background 
spectra were extracted from regions as described in Sect. 2.1. 
We limited the energy band of the pn camera to 0.3 - 2.2 keV, 
while the MOS cameras were analysed between 0.15keV and 
2.2 keV. 

To increase the S /N and use standards-minimization tech- 
niques in XSPEC 12.7.0, spectra were rebinned using the SAS 
task specgroup. The spectral bins were grouped according to a 
minimum S /N of three and we took care to oversample the in- 
strument energy resolution at a given bin by less than a factor of 
three. Our preliminary analysis had shown that the inclusion of 
data from the MOS2 camera considerably worsens the spectral 
fits and biases the parameter determination; we therefore only 
considered the pn and MOS 1 cameras, which were fitted simul- 
taneously to more tightly constrain the spectral parameters. 

To account for interstellar absorption, we adopted the pho- 
toelectric absorption model and elemental abundances of Wilms 
et al. (2000, tbabs in XSPEC). We note that the source spectral 
parameters in Paper I were derived by adopting a different abun- 
dance table (that of Anders & Grevesse 1989), which produces a 
~ 15% higher value of the equivalent hydrogen column-density 
and a ~ 10% softer blackbody temperature than those derived 
with Wilms et al. abundances (c.f. note 3 ). 

The results of our spectral fits are given in Table 4. We 
list for each fit (numbered [1-14] as reference for the text) the 
reduced chi-squared xl, degrees of freedom and correspond- 
ing null-hypothesis probability, the equivalent hydrogen column 
density Nu, model-dependent parameters, and the unabsorbed 
source flux Fx in the 0.1 - 12keV energy band. Unless otherwise 
noted, the errors in the parameters represent the 1 cr confidence 
levels. 

We tested several models and found that the emission of 
J 1046 is consistent with previous analysis based on archival 
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Model Xv (d.o.f.) Null Hyp. Parameter ~F\ Comment 







(%) 


(10 21 cm" 2 ) 








(erg s 1 cm 2 ) 




[1] bbody 


1.66(61) 


0.091 


2 61 +U 22 


kToc 


138(3) 


eV 


6.8!^ x 10-" 


Wilms et al. abundances 


[2] bbody 


2.22(61) 


1.3 x 10~ 5 


2.76(23) 




152(3) 


eV 


5.9(4) x 10 


Hamaguchi et al. abundances 


[3] bbody 


1.21(60) 


13 


q q i +0.28 
-0.27 




124(4) 


eV 


1 f >A/n\ w 1 a— 12 
1.20(22) X 10 


Zo = 1.66(12) (solar) 


[4] nsa B 


1.64(61) 


0.114 


^" iJ -0.16 




4 - 59 -ai5 x 1(P 


K 


3.3 + _\ l x 10-' 2 


M* = 1.4M ,R* = 10km 


[5] nsa Bo 


1.70(59) 


0.067 


4.25(18) 




4.5^1 x 10 5 


K 


33 + _\l x 10- 12 


M = 1.3+J5 3 M ,R ~ 10km" 


[6] nsa B 12 


1.80(61) 


0.013 


4 20 +0 - 22 




7.1(3) x 10 s 


K 


2.8^ x 10- 12 


M* = 1.4M ,R* = 10 km 


[/ J nsd d 12 


i .^rjyoy } 


1.673 


? 1 5 +0 ' 27 
^ -0.20 




3 n +0 - 7 y i n 6 

5 




4+5 v in-' 3 


M ~ 9 M fi* ~ 7 hn ^ 
JV1 ~ Z IVIq, i\ ~ / Kill 


[8] nsa B 13 


1.98(61) 


7.6 x 10 


A ££+0 21 

4.66_ 04 




£ o i +0 27 w 1 r\5 
O.B 1 Iq 3 X 10 


K 


3.9! 3 6 x 10-' 2 


M* = 1.4 M , R* = 10km 


[9] powerlaw 


2.37(61) 


9.1 x 10"' 


9.9(4) 


r 


8.58(26) 




8(5) x 10"' 


Wilms et al. abundances 


[10] bbody-gauss 


1.05(57) 


38 


2.7(3) 


kToo 


1351 


eV 


7.5^ x 10-' 3 


Z = 1.64(15) (solar) 










£ 


1 36 +003 


keV 


o- = 0.14(lo-),EW = 91eV 


[11] bbody-2*gauss 


1.01(58) 


46 


2.61* 


kT m 
ei 

<?2 


134.9(2.4) 

o 6ir 0020 

u.uii _ Q018 

1.35(4) 


eV 

keV 

keV 


7.3(4) x 10- 13 


Wilms et al. abundances 
cr = 0.1*,EW = 71eV 
<t = 0.1*, EW = 86eV 


[12] bbody-2*gauss 


1.00(57) 


48 


2.2(3) 


kToo 


141(6) 


eV 


5.7+> 2 x 10- 13 


Wilms et al. abundances 










ei 


597+ - 023 


keV 


o- = 0.1*,EW = 71eV 










£2 


1.35(3) 


keV 




o- = 0.1*,EW= 104eV 


[13] bbody+bbody 


1.17(59) 


17 


5.2(6) 




37(4) 


eV 


1.3! 27 8 x 10- 10 


Rf ~ 640 km c 










*r 2 » 


119(4) 


eV 


K~ ~ 6km r 


[14] bbody+pow 


1.66(61) 


0.091 


2.61* 


r 


137.9!|- 2 
1.7* -2.1* 


eV 


6.8(3) x 10- 13 


Wilms et al. abundances 

Ff = (6.8- 7.3) X 10- 15 (3o-) d 



Notes. Parameters marked with a star are held fixed during fitting. 

(fl) Stellar radius unconstrained by the model. (b) Stellar mass and radius unconstrained by the model. (c) Computed for distance d = 2.3 kpc. 
{d) Upper limit to the unabsorbed flux of the power-law component, in units of erg s~' crrr 2 (0. 1 - 12 keV). 



data 3 . However, the much higher quality of our statistics revealed 
that a single component is hardly satisfactory in describing the 
source energy distribution. The best-fit single-component model 
consists of either a blackbody or a non-magnetized neutron star 
atmosphere. The null-hypothesis probability, which gives the 
likelihood of having a larger value of xl by chance, assuming 
that the model is correct, is low - at best lower than 2% for all 
single-component models ([1-9] in Table 4). The best-fit power- 
law model [9] is significantly worse than the thermal models, 
with a null-hypothesis probability <k 1%. The photon index 
T ~ 9 is unreasonably steep and the model is four times more 
absorbed than the thermal ones. Other models - bremsstrahlung 
and optically thin thermal plasma (Raymond & Smith 1977) - 
show even poorer fits. We discuss below the results of the ther- 
mal models in more detail. 

The best-fit absorbed blackbody [1] has kT^ = 138(3) eV 
and /V H = 2.61+° ^ x 10 21 cnr 2 , with xl ~ 1-7 for 61 degrees 
of freedom. The radiation radius, as seen by an observer at infin- 
ity, is R m ~ (3 - 5)(if/(fcar) 1 km, where the distance to J1046 is 
normalised to that of 77 Car, dcm = 2.3 kpc (Smith 2006b). In this 
estimate, we included the flux uncertainty from the abundance 
tables. Although this radius is smaller than that of a canonical 
neutron star, it agrees with those measured for the M7: their red- 
shifted radiation radii, as derived from X-ray blackbody fits and 
distance estimates, are in the range of, roughly, 2 km to 7 km 
(using the distance estimates and source parameters from e.g. 
Kaplan 2008; see also Sect. 5.3). 



3 The best-fit blackbody model gives kT m = 122(4) eV and A^h = 
3.20+q 27 x 10 21 cm -2 , with the abundance table of Anders & Grevesse 
(1989). 



A neutron star atmosphere model (Zavlin et al. 1996; Pavlov 
et al. 1995, nsa) describes the data statistically nearly as well [4- 
8]. A stellar mass of M ns = 1.4M and radius of R em = 10 km 
were first assumed [4,6,8], then allowed to vary to check for 
an improved fit [5,7]. The magnetic field strength was held 
fixed during fitting, and we tested the values B — 0G [4-5], 
B = 10 12 G [6-7], and B = 10 13 G [8]. The best-fit models typi- 
cally fit a 60% higher /V H than the blackbody and a much softer 
temperature, kTco ~ 30 eV. For the typical absorption in the di- 
rection of the Carina Nebula, the implied luminosity distance is 
implausibly small, d ~ 220 pc. Poorer fits were found when the 
magnetic field was included and with increasing field intensity. 
Although the fit was generally not improved when the neutron 
star mass and radius were allowed to vary, we found a statisti- 
cally tighter fit for a magnetic field strength of B — 10 12 G and 
with freely varying M ns ~ 2 M and R m ~ 7 km (null-hypothesis 
probability of 1.7%; [7]). In this case, A^h and kT^ were more 
similar to the most closely fitting blackbody values. However, 
the stellar mass and radius were largely unconstrained by the 
model, and the implied source distance is much larger than ex- 
pected from the measured /V H (d ~ 7.5 kpc). 

The best-fit blackbody/neutron star models invariably show 
residuals around 0.6 - 0.7 keV and 1.3 - 1.4keV (see Fig. 5; 
left). We note that these features are independent of the choice 
of source and background regions used to create the spectra (thus 
exclude improper background subtraction), of the EPIC instru- 
ment (they are seen in the data from the three EPIC cameras, 
but are more pronounced in pn), as well as of the choice of 
cross-sections and abundance tables available in XSPEC. Before 
adding yet more complexity to the spectral model, we investi- 
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Energy (keV) Energy (keV) 

Fig. 5. Results of spectral fitting of source J1046. We show the pn data and folded best-fit blackbody model, with residuals. Left. Simple blackbody. 
Right. Best-fit with oxygen overabundance and a Gaussian absorption line at energy e ~ 1.35 keV. 



gated whether the residuals could be due to a local elemental 
overabundance in the Carina Nebula. 

Recent studies (e.g. Townsley et al. 2011a, and references 
therein) have investigated the properties (morphology, emissiv- 
ity, and elemental abundances) of the hot plasma known to per- 
meate the Carina Nebula. The interaction of the plasma with cold 
molecular gas is likely responsible for the diffuse soft X-ray 
emission first detected with the Einstein Observatory (Seward 
et al. 1979; see Sect. 5.2). Its composition and, in particular, en- 
hanced abundances of key elements indicate whether this soft 
emission originates in either the wind of the current population 
of more than 70 massive stars (Smith 2006a) and/or in one or 
more supernova explosions (Hamaguchi et al. 2007). 

We tested the hypothesis of a local elemental overabundance 
using the variable photoelectric absorption model vphabs con- 
volved with a blackbody in XSPEC. We first attempted to fit the 
spectra by fixing the elemental abundances at the best-fit values 
found by Hamaguchi et al. (2007), with freely varying Nu and 
kT [2]. However, the result was considerably worse than that 
obtained with solar abundances from Wilms et al. Moreover, the 
residuals remained roughly at the same energies. We then tested 
for the overabundance of all elements with significant transitions 
in the range of 0.2 keV to 2keV (C, N, O, Ne, Mg, . . . ). We re- 
stricted our test to abundance values in the range 1 < Z < 5 
in solar units to avoid both subsolar and arbitrarily high metal 
abundances (see e.g. Townsley et al. 201 la, for a discussion). 

Interestingly, we found that, whereas most of the elemen- 
tal abundances (e.g. C, N, S) remained insensitive to the fit and 
others (most notably, Mg) were pegged to the hard end of the 
allowed interval in abundance, the only element that was tightly 
constrained was oxygen, with Zq = 1.66(12) in solar units [3]. 
By allowing Zq to be extrasolar, we obtained a much closer 
agreement between data and model for energies below 1 keV. 
However, an inspection of the residuals revealed that the feature 
at energy ~ 1 .3 keV remained in spite of the improved fit and 
could only be accounted for with a Gaussian line in absorption. 
We indeed verified that elements with lines around this energy 
tended to be arbitrarily overabundant in the model. 

The inclusion of a Gaussian line in absorption improves the 
fit considerably [10]. Allowing all parameters to vary freely, 
we obtained a good fit with xl ~ 1 an d a much higher null- 
hypothesis probability. The best-fit model and its residuals can 



be seen in the right panel of Fig. 5. This model - with Nu and kT 
consistent with the (simple) blackbody values, an oxygen over- 
abundance of Zo = 1.64(15), and a line in absorption at en- 
ergy e = 1.36^QQ5keV (<x = 0.14keV and equivalent width of 
EW = 91 eV) - reproduces the observed spectral energy distri- 
bution of J1046 very well. We verified that the results (e, Zq) are 
insensitive to any changes in A^h- 

For completeness, we tested other multicomponent models, 
namely a double blackbody [13], a blackbody with multiple 
Gaussian absorption lines (with fixed solar abundances, [11- 
12]), and a blackbody with a power-law tail [14]. The fit results 
are also listed in Table 4. All parameters in the double black- 
body [13] were allowed to vary during fitting. The softer black- 
body temperature kT\ was restricted to between 5 eV and 120 eV 
and the hotter one kT2 to between 0.1 keV and 1 keV. The best 
fit (xl ~ 1.2) gives temperatures of 37 eV and 119eV. Relative 
to the single-blackbody fit, the addition of the extra component 
provides a closer fit to the softer part of the spectrum, elimi- 
nating the "excess" detected previously and lowering the value 
of x 2 - However, the two features at energies ~ 0.6 keV and 
~ 1.3keV discussed before remain. Given the substantial inter- 
stellar medium absorption of the model (Nu ~ 5 x 1 21 cm -2 ), the 
normalisation of the soft component is forcibly very high, two 
orders of magnitude higher than the hot component. The red- 
shifted emission radius, at a distance of 2.3 kpc, is R™ ~ 640 km 
and = 6 km. This argues against the physical viability of the 
model. 

To test for the presence of multiple absorption lines, we 
adopted the following strategy: first, we fitted the data hold- 
ing the column density fixed at the best single-blackbody value, 
Nn = 2.61 x 10 21 crrT 2 , also keeping the Gaussian cr parameter 
fixed at a value of 0. 1 keV [11]. We then allowed to vary either <x 
or A^h- Finally, we fitted the data allowing all parameters (but the 
equivalent widths) to vary freely [12]. For all fits, the blackbody 
temperature was restricted to within lOeV and 300 eV and the 
line energy to within 0.1 keV and 2keV. Spectral fits show good 
results with two absorption lines invariably at energies around 
0.58 - 0.61 keV and 1.35 - 1.38 keV. The presence of the lines 
does not lead to any significant changes in the best-fit values of 
kT and Nu- When Nu varies freely, the blackbody is hotter rela- 
tive to the single-component temperature (kT ~ 140 eV) and the 
column density is lower, at Nu ~ 2.2 x 10 21 crrT 2 [12]. 
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Table 5. Description of the optical ESO-VLT observations of J1046 



Night 


Period 


Exposures 




FWHM 


Airmass 








(s) 


(arcsec) 




2009-02-21 


P82 


3 x V 


3300 


0.50 


1.22 


2010-04-07 


P85 


6 x Ha 


6264 


1.16 


1.23 


2010-06-06 


P85 


3 xHa 


3132 


0.67 


1.26 



Notes. The exposure times, seeing, and airmasses are averages per filter. 
(a) Total exposure time per filter. 



We note that Hamaguchi et al. (2009) reported evidence of 
absorption features in the combined spectra of J 1046, namely a 
slight dip around 0.9 keV and a strong one at 0.6 keV. As noted 
by the authors, the latter is likely to be produced by edge ab- 
sorption of oxygen in the interstellar medium or in the detector 
response. We found no evidence in A09 data to confirm the dip 
at 0.9 keV, which might have been produced during the combi- 
nation of both the spectra and calibration files of the three EPIC 
detectors taken at various off-axis angles and with different in- 
strument configurations. 

The blackbody model in combination with a power law [14] 
was tested to derive the upper limits to the detection of a non- 
thermal component extending towards higher energies, as usu- 
ally seen in the emission of middle-aged pulsars dominated by 
soft thermal components (age ~ few 10 5 yr, e.g. the "Three 
Musketeers"); in these objects, power-law tails with Y — 1 .7-2. 1 
were detected at very low flux levels (0.3 - 1.7% of the source 
luminosity; De Luca et al. 2005). To derive upper limits to the 
detection of a hard non-thermal component, we fitted a power 
law with photon indices r = 1.7 and F = 2.1, in addition to the 
dominant blackbody component. The photon index and a col- 
umn density of Nn = 2.6 x 10 21 ctrT 2 were kept fixed during the 
fitting. We found that the best-fit blackbody parameters were un- 
affected by the inclusion of the power-law component. Relative 
to the single-blackbody model, the fit did not improve and the 
non-thermal component was found to contribute at most 1-1.1% 
(3cr confidence level, 0.1 - 12keV range) of the unabsorbed flux 
of the source. 

3. ESO-VLT observations 

Optical follow-up observations of J 1046 with ESO-VLT and 
SOAR in 2007 and 2008 revealed no counterparts brighter than 
mR = 25, mB = 26, and my = 25.5 (2.5<x; Paper I). To derive a 
deeper limit to the brightness of the optical counterpart, we ob- 
tained further ESO-VLT V-band exposures in 2009. The V band 
is the optimal choice for deep imaging because it excludes some 
of the strongest lines found in H n regions (Balmer emission as 
well as [N n] and [S n] emission lines), thus minimizing con- 
tamination from bright nebular lines. Additionally, we obtained 
deep Ha imaging with ESO-VLT in 2010, in order to investigate 
a possible B aimer-dominated nebula near the position of the X- 
ray source. 

3.1. Observations and data reduction 

The VLT 2009 and 2010 observations were obtained under pho- 
tometric and dark sky conditions for a total ~ 3.53 h of ob- 
serving time (see the log of observations in Table 5). The to- 
tality of the proposed observing time was executed in two dif- 
ferent observing periods (P82 and P85). On both occasions, 
one of the FOcal Reducer low/dispersion Spectrographs (FORS; 
Appenzeller et al. 1998) was used. 



The FORS detector consists of a mosaic of two 2kx4k CCDs 
of pixel size 15yum. In the 2010 observations (P85), we used 
FORS2, the camera mounted on UT1, which is optimized for 
the red band (MIT CCDs), while the 2009 observations of P82 
were carried out with FORS1 (E2V CCDs) at UT2. The spectro- 
graph FORS1 is most sensitive to the blue range. The detector 
provides imaging on a pixel scale of 0.25" pixeL 1 (field-of-view 
of 6.8' x 6.8') using standard-resolution collimator and default 
readout mode (2x2 binning). Dithering patterns with offsets of 
5" were chosen for the observations (with 3 x V and 9 xHa expo- 
sures; see Table 5). The seeing was excellent during the 2009 V 
observations; unfortunately, it was above the constraint of 0.8" 
FWHM during the observing night of April 2010, when most 
(67%) of our P85 observations were executed. 

Data were reduced in two steps. Firstly, we applied 
IRAF v2.14 (Tody 1986) to process raw frames; secondly, 
EsoRex 4.3.5 4 and ESO-VLT recipes were used in order to recre- 
ate the instrument pipeline. For both sets of software, stan- 
dard procedures were adopted and carried out by treating each 
FORS CCD chip independently. The flat-fielded scientific expo- 
sures were combined in order to remove both cosmic rays and 
bad/hot pixels and increase the S /N of the data. We verified that 
both (IRAF/EsoRex) reductions yielded very similar results in the 
analysis, which was also true when compared to the results ob- 
tained with the pipeline-reduced scientific frames distributed by 
ESO-VLT. 

3.2. Data analysis 

The transformation of the my instrumental magnitudes to the 
standard photometric system was done using the zero-point, ex- 
tinction coefficient and colour term obtained on the night of 
the observation, as provided by the ESO-VLT calibration web- 
pages 5 . The Ha images were flux-calibrated using observations 
of a standard star observed in the same night of the observa- 
tion (LTT4816; Hamuy et al. 1994). The astro me trie calibration 
of the stacked science frames was performed using the USNO- 
B1.0, 2MASS, and GSC2 catalogues and the GAIA 4.2-1 soft- 
ware 6 . In general, our astrometric errors are of ~ 0.15" or better. 
We adopted standard PSF fitting, as implemented in the daophot 
(Stetson 1987) package for the IRAF environment, to measure 
magnitudes. 

The deep V observation failed to detect the optical counter- 
part to J1046: as can be seen in the optical image of Fig. 6, no 
object lies within ~4.3" (>5<x) of the source position. To test 
the detection limit specifically in the region covered by the error 
circle of the X-ray source, we added images of synthetic stars of 
progressively fainter magnitudes to a small section of the com- 
bined V image and analysed the combined images using PSF 
fitting (see Pires et al. 2009a, for details). We defined the limit- 
ing magnitude of the 2009 observations as my = 27.6(5) (2cr), 
which corresponds to the detected magnitude of the faintest syn- 
thetic star successfully measured in the image. Taking the lim- 
iting magnitude as an upper limit to the brightness of the opti- 
cal counterpart of J 1046, the implied X-ray-to-optical flux ra- 
tio is log(F x /^v) > 3.8, corrected for photoelectric absorption 
and interstellar extinction. The de-reddened V flux is derived us- 
ing the X-ray absorption (Sect. 2.4) and the optical extinction, 
Ay ~ 1.51, computed with the Predehl & Schmitt (1995) rela- 
tion. 

4 http://www.eso.org/sci/software/cpl/esorex.html 

5 http ://w w w. eso . org/ observing/df o/quality 

6 http://star-www.dur.ac.uk/~pdraper/gaia/gaia.html 
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Fig. 6. V (left) and Ha (right) images of the field of source J1046. An 
confidence level error circle on the position of the source is of size 0.86". 

3.3. Searching for a bow-shock or photoionisation nebula 

The dense environment of the Carina Nebula and the possibly 
lowly-ionised local medium surrounding J1046 favour the for- 
mation of a bow-shock or of a photoionised nebula that could be 
detected in the optical. The fraction of neutral hydrogen in the 
nebula depends on the distance from the sources of ionising radi- 
ation, mostly OB stars present in the young stellar clusters Tr 16 
(where 77 Car is) and Trl4, and on the dumpiness of the local 
medium. J 1046 is located in a region of low diffuse background 
relative to other parts of the nebula, and is apparently behind an 
ionisation front from Tr 16 (see Fig. 7). Therefore, most of the 
local interstellar medium may still be neutral, which favours the 
detection of an Ha signature around the position of the presumed 
rapidly-moving neutron star. 

A Balmer-dominated nebula was seen around a number of 
radio pulsars (e.g. Gaensler & Slane 2006) and one (possibly 
two, see Motch et al. 2005 and the case of RX J1605. 3+3249) 
of the M7, RX J1856.5-3754 (van Kerkwijk & Kulkami 2001). 
Two possible scenarios are invoked to explain the presence of 
such a nebula: bow-shocks (e.g. Chatterjee & Cordes 2002) 
and photoionisation by X-rays (Blaes et al. 1995). The shape 
and emissivity of the nebula can provide invaluable information 
about the INS velocity, direction of motion, and rotational en- 
ergy losses, allowing one to constrain its age and evolutionary 
state. 

Motivated by the above, we estimated the expected Ha flux 
from J1046 for both scenarios, considering angular (projected) 
separations of the apex of the nebula relatively to the neutron 
star position (stand-off radius Rq) in the range from 0.5" to 5". 
Simulations using the photoionisation code XSTAR 7 predicted 
an Ha flux (taking into account the interstellar extinction) of 
fua ~ 10~ 17 - 10~ 16 erg s _1 ctrT 2 arcsec~ 2 , for a local density of 



7 http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 
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inverted colour map is used, i.e. brighter objects are darker. The 90% 
The images were smoothed using a Gaussian filter of one pixel in size. 




Fig. 7. ESO-MAMA wide-field R band image of the area around J1046 
(marked by a cross), in the direction of the Carina Nebula. The young 
OB associations Tr 16 and Tr 14, as well as 77 Car, are located NW of the 
source, towards the bright diffuse emission. 



50 cm (e.g. Paper I) assuming the spectral properties of J1046 
(Sect. 2.4). The resulting flux range is consistent with the ex- 
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pectations of Blaes et al. (1995) for a similar neutron star lumi- 
nosity and medium density, and with that measured for the M7 
RXJ1856.5-3754. 

Alternatively, in the bow-shock scenario, a neutron star of 
velocity v ns = 10 - 100km s~' moving through a uniform 
medium of density n = 50 cm 3 , where (conservatively) half of 
the hydrogen atoms is still neutral, requires a spin-down lumi- 
nosity of E > 10 34 ergs _1 to create bow-shocks with stand-off 
radii in the range Rq = 0.5" - 5" (see the formulation of e.g. 
Cordes et al. 1993). These shocks correspond to an Ha flux 
of fu a ~ 10~ 15 - 10~ 14 ergs~' cirT 2 arcsec~ 2 , accounting for 
absorption. If, conversely, the spin-down luminosity of J 1046 
is lower, more typical of those of the M7 (E ~ 10 31 ergs -1 ), 
the Ha nebula is then expected to be less bright, fu a ~ 4 x 
10~ 17 ergs~' cirT 2 arcsec~ 2 , and (at d = 2.3 kpc) nearly angu- 
larly coincident with the position of the neutron star. In any case, 
these figures indicate that the flux in the Balmer line may be high 
enough to allow detection with ~ 3 h of VLT observations, for 
both (photoionisation/bow-shock) scenarios. 

Previous Ha imaging of the field of the X-ray source, 
obtained with SOAR in 2008, was badly affected by fring- 
ing (Paper I). Our new VLT observations do reveal some ex- 
tended structure in the vicinity (within 1") of the position of 
the X-ray source (Fig. 6; right). The Ha flux is f Ha ~ 3.6 X 
10~ 15 ergs~' cirT 2 arcsec~ 2 - at least one order of magnitude 
brighter than that expected if the nebular emission was created 
by X-ray photoionisation as previously discussed, but still for- 
mally consistent with the expectations of the bow-shock scenario 
forZ? > 10 34 ergs~' (typical of young rotation-powered pulsars). 
However, given the presence of similar diffuse emission in other 
parts of the nebula, it is unlikely that the observed structure is 
indeed related to the neutron star. 

4. Gamma-ray analysis 

Our measurement of the neutron star spin period permits us to 
search for pulsed emission also in y-rays, using data from the 
AGILE and Fermi-LAT missions. This is particularly interest- 
ing given the large time spans of data that can be used to con- 
strain the period derivative of the pulsar. Although no plausi- 
ble y-ray counterpart is present within > 5cr of the position of 
the X-ray source in either the AGILE Bright Source Catalogue 
(Pittori et al. 2009) or the Second Fermi-LAT Catalogue (2FGL; 
The Fermi-LAT Collaboration 201 1), a search for a periodic sig- 
nal can be more sensitive than spatial detections (see Zane et al. 
201 1, and references therein, for a discussion). 

Data reduction was performed by making use of the software 
specific to each telescope and the whole available data span of 
both satellites. Public data from the AGILE mission were ac- 
cessed via the Agenzia Spazionale Italiana Science Data Center 
(ASDC) Multi-Mission Interactive Archive 8 . The analysis cov- 
ers the time span from 2007 July 13 (MJD 54294) to 201 1 April 
5 (MJD 55656), which amounts to 169,289 photons extracted 
around 5° (in order to match the point spread function of the in- 
struments at energies E > 100 Me V) from the position of J 1046, 
as reported in Table 2. We applied the AGILE pulsar pipeline 
to perform data selection and barycentering. For details of the 
AGILE observing strategy, timing calibration, and y-ray pulsar 
analysis, we refer the reader to Pellizzoni et al. (2009a,b). 

The public archive of the Fermi Science Support Center 9 was 
accessed to retrieve Fermi-LAT data. Again defining a region- 

8 http://www.asdc. asi.it/mmia/index.php ?mission=agilemmia 

9 http://fermi.gsfc.nasa.gov/cgi-bin/ssc/LAT/LATDataQuery.cgi 



of-interest of 5° around the position of the source, 760, 327 pho- 
tons with energy higher than 100 Me V and registered between 
2008 August 5 (MJD 54683) and 201 1 March 30 (MJD 55650) 
were extracted in the timing analysis, using the Fermi-LAT 
Science Tool 10 gtselect. Good-time intervals and data filters rec- 
ommended by the Fermi-LAT Team were also applied with gtmk- 
time. Finally, the photons were barycentred using the gtbary tool. 

We searched for periodicities using as a starting point the 
spin period of the source observed in X-rays (Sect. 2.3). The two 
data sets (AGILE/Fermi-LAT) were used separately and simul- 
taneously. A pulsed signature was searched in frequency domain 
close to the nominal frequency v, = 53.645804 Hz and around 
it, using a frequency step of the order 1/T S (where T s is the time 
span of the analysed data set). To look for significant pulsations, 
we applied a Pearson x 2 test to the ten-bin pulsed profile result- 
ing from each search trial. We also applied a (bin-independent) 
Z 2 test (Buccheri et al. 1983) to each analysed data set. 

It is worth noting that the restricted information available for 
J1046 greatly enhances the search in terms of parameter space. 
The number of trials, given the long T s for both the AGILE and 
Fermi-LAT data sets, significantly increases, even when search- 
ing only around 1 cr in frequency from the nominal value. As a 
first approach, the first and second period derivatives were set 
to zero and not allowed to vary. A wider search was then per- 
formed across a small range of period derivatives. Furthermore, 
to extend the search to a broader frequency range (around 10<x 
of v,), the time span of the analysis was restricted to cover only 
the photons registered in December 2010, when the X-ray obser- 
vations took place. 

No significant periodicity (> Act) was found. Although tenta- 
tive pulsations result from individual searches, overall inconsis- 
tencies prevent a secure detection from being claimed in y-rays. 
Hence, we set a Act upper limit of F y < A x 10~ 8 ph s" 1 cirT 2 at 
energies above 100 MeV, for a pulsed y-ray signal from J1046. 

5. Discussion 

We first summarize our results and then discuss the origin of 
J 1046 in relation to the former generation of massive stars in 
the Carina Nebula and its possible evolutionary state. We in- 
terpret our findings in the light of the observed properties of 
the currently known neutron star population and compare J 1046 
to other peculiar objects. The intriguing possibility of detect- 
ing J 1046 with the current and future generation of gravitational 
wave detectors is briefly discussed. Finally, we consider the pos- 
sible number of Galactic neutron stars and the prospects of find- 
ing new thermally emitting sources with the eROSITA mission 
(Predehl et al. 2010), planned to be launched in 2013. 

5.1. Summary of results 

We have presented the results of the first dedicated observational 
campaign designed to investigate the properties of the isolated 
neutron star 2XMM J 104608. 7-594306. Our analysis is based on 
new observations obtained with the XMM-Newton Observatory 
and the ESO Very Large Telescope, as well as publicly avail- 
able y-ray data from the Fermi Space Telescope and the AGILE 
Mission. Our observational campaign has confirmed expecta- 
tions based on a previous analysis of archival data and revealed 
a unique and peculiar object. 

The new XMM-Newton observation alone collected a fac- 
tor of two more source counts than all previous X-ray obser- 

10 Available at http://fermi.gsfc.nasa.gov/ssc/data/analysis/software 
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vations (that serendipitously detected the neutron star in the past 
12 years) did collectively. The higher quality statistics and higher 
time resolution of the pn camera in small window mode, has per- 
mitted the measurement of the period of the very likely neutron 
star spin, at an unexpectedly rapid rotation of P ~ 18.6 ms. The 
detection of the periodicity, corresponding to a pulsed fraction 
of pi ~ 14%, has a chance of 1 in 16,000 of being spurious. 
No other periodicity is significantly detected (> lcr) across a 
very broad period range, P = 0.01 14 - 10000 s in pn data. The 
analysis of the three EPIC cameras together more tightly con- 
strained the periodicities with pf > 10% (3cr) and P > 0.6 s 
(0.3 - 2.5 keV), thus considerably improving previous limits in 
this frequency range. 

The XMM-Newton observation also confirmed the purely 
thermal nature of the source. The energy distribution is most ac- 
curately decribed by a simple blackbody of temperature kT m ~ 
135eV, with significant deviations (absorption features) around 
energies 0.6 keV and 1.35 keV. While the first feature is more 
likely related to a local oxygen overabundance in the Carina 
Nebula or along the line-of-sight, the second one can only be 
accounted for by an additional spectral component, modelled as 
a Gaussian line in absorption. The source emission is very soft 
with practically no counts above 2 keV; in contrast to standard 
rotation-powered pulsars and magnetars, a non-thermal compo- 
nent extending towards higher X-ray energies is detected in- 
significantly (< 3cr) at flux levels higher than 1% of the source 
luminosity (L x = 4.7+^ x 10 32 ergs _1 , 0.1 - 12keV). Similarly, 
we found no evidence of a two-temperature spectral model, as is 
often the case for thermally emitting radio-quiet neutron stars lo- 
cated near the centre of supernova remnants (CCOs; see de Luca 
2008, for a review), some transient AXPs, or Calvera (Sect. 5 .3). 

The new ESO-VLT observations set an even deeper limit 
(relative to previous optical follow-up) to the brightness of the 
optical counterpart of J1046, m v > 27 (2cr); the optical flux is 
therefore at least 6300 times fainter than the X-ray emission of 
the neutron star, confirming the compact nature of the source. 
Although a diffuse Ha excess is present near the position of the 
X-ray source, similar background emission in other parts of the 
nebula indicates that it is unlikely to be related to the neutron 
star. 

Our analysis of public y-ray data has revealed no signifi- 
cant detection. Given the still unknown and unconstrained pe- 
riod derivative of the pulsar (see our discussion below), it is dif- 
ficult to say whether any should be expected. By analogy with 
other INSs, it is likely that the pulsar spin-down luminosity is 
too small, given its estimated distance, to observe y-ray emis- 
sion with the present sensitivity. 

5.2. Former generation of massive stars in the Carina Nebula 

The Great Nebula in Carina is an ensemble of young and rich 
stellar clusters, highly structured dust clouds and pillars, cavi- 
ties, and bubbles created by the dramatic interaction of ionisation 
fronts with scattered molecular material (see Smith & Brooks 
2008; Townsley et al. 201 lb, for reviews). Star formation is cur- 
rently underway, as revealed by the presence of e.g. pre-main- 
sequence objects, microjets as well as many candidate young 
stars detected in the infrared (Smith et al. 2010a,b). Moreover, 
populations of stars with ages of from 5 to lOMyr, as in Tr 15 
(Wang et al. 201 1), support evidence that star formation has been 
ongoing for the past several million years - i.e., a sufficiently 
long time for the most massive objects belonging to the for- 
mer stellar generation in the complex to have already ended their 
lives in core-collapse events. 



Although for instance Smith (2006a) and Smith & Brooks 
(2008) argue that past supernovae are not needed on energetic 
grounds alone to create the observed cavities and superbubbles, 
there is growing evidence that at least one explosion has oc- 
curred in Carina: (z) the detection of an unusual and broad emis- 
sion feature at 22 pm, related with recently synthesized dust, 
which has only been observed in the Cassiopeia A supernova 
remnant and in starburst galaxies (Chan & Onaka 2000); (if) en- 
hanced abundances of iron and silicon in the southern part of the 
nebula relative to other regions, which indicate that the metal en- 
richment originated in a supernova rather than winds of evolved 
stars (Hamaguchi et al. 2007; Ezoe et al. 2009); {in) the presence 
of high-velocity expanding structures, which are predominantly 
located in front of main sources of ionisation, such as the stellar 
cluster Tr 16 (Walborn et al. 2007); (z'v) the evidence of charge 
exchange as the origin of Carina's diffuse soft X-ray emission, 
since a rarefied hot plasma resulting from an old supernova is 
a necessary condition for generating such a process (Townsley 
etal. 2011a). 

The discovery of J 1046 provides further support for the past 
supernova hypothesis (Hamaguchi et al. 2009; Townsley et al. 
201 la). The column density derived from the X-ray spectral fits 
is indeed consistent with the source being located in the Nebula; 
its softness excludes an extragalactic object, given the Galactic 
absorption along the line-of-sight (7Y H = 1.35 x 10 22 cirT 2 ; 
Dickey & Lockman 1990). Townsley et al. (2011a) suggested 
that J 1046 might have been the first discovered neutron star and 
that others are still hidden in the Nebula, proposing an additional 
six X-ray sources without any evident optical or infrared coun- 
terparts as possible INS candidates (see their Table 2). In prac- 
tice, however, the faintness of these candidates (observed fluxes 
of fx ~ 10~ 15 - 10~ 16 ergs~' crrT 2 in the 0.5 - 8keV energy 
band; Broos et al. 2011) requires unrealistically deep follow-up 
optical observations in order to confirm their tentative compact 
nature: observed upper limits in range m B = 29-30 are needed to 
reach an X-ray-to-optical flux ratio of 100, which is still not high 
enough to discard either cataclysmic variables or background ac- 
tive galactic nuclei. 

If J 1046 is indeed a remnant of a former population of mas- 
sive stars in Carina, the source is still somewhat young; younger 
than several 10 6 yr, considering typical main-sequence phases of 
0.5-5Myr. 



5.3. Possible nature and relation to other neutron stars 

The new results of our observing campaign targeting J 1046 can 
be used to understand how peculiar groups of isolated neutron 
stars relate to each other, as well as to the bulk of the nor- 
mal radio pulsar population. We therefore compare the spec- 
tral and timing properties of possible "missing links" - J 1046, 
RRAT J1819-1458, and Calvera - with those of other ther- 
mally emitting INSs, namely the M7, magnetars, CCOs, and 
millisecond pulsars (MSPs). We first briefly review the prop- 
erties of CCOs and MSPs in the next subsection (for magne- 
tars, the M7, Calvera and RRATs, we refer to Sect. 1, and ref- 
erences therein). In the discussion, we assume that the pulsar 
spin down is due to magnetic dipole braking, i.e. dipolar mag- 
netic field fidip = 3.2 x 10 19 (PP) 1/2 G, spin-down luminosity 
E = 4.5 x 10 46 (Pf- 3 )ergs- 1 , and r ch = P/(2P)s, where r ch 
is the pulsar characteristic time. 
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Fig. 8. Blackbody luminosity versus emission radius (left) and emission radius vs. neutron star spin period (right) for several populations of 
thermally emitting INSs (see text). 



5.3.1 . Properties of thermally emitting INSs 

Millisecond pulsars are old and low-magnetized (r c h = 0.1 - 
10 Gyr and Z?di P ^ 10 10 G) neutron stars, mostly observed in both 
binary systems and globular clusters. Accretion in a binary sys- 
tem is thought to be responsible for "recycling" these old neutron 
stars, first by spinning up the pulsar to millisecond periods so ra- 
dio emission is turned on again, second by heating areas of the 
surface - polar caps or hot spots - which can then emit X-rays. 
The bulk of the X-ray radiation is usually thermal and originates 
from the hot spots. Power-law tails, or purely non-thermal X-ray 
emission, are also often observed in these objects. Interestingly, 
Fermi-LAT has also revealed that MSPs are a major contributor 
to the Galactic y-ray source population, providing insights into 
the magnetospheric emission processes (e.g. Abdo et al. 2009, 
2010). 

Of the seven confirmed CCOs, at least three fit the "anti- 
magnetar" scenario (Halpern & Gotthelf 2010; Gotthelf & 
Halpern 2010), in which these objects are young (ages 10 3 - 
10 4 yr), low-magnetized (B < few 10 n G) neutron stars, born 
rotating at a period similar to its present value (of hundreds of 
milliseconds). Spin-down values and upper limits are very low, 
confining anti-magnetars to a "transition" region of the P - P di- 
agram, in-between those occupied by normal rotation-powered 
pulsars and MSPs. The energy distribution of CCOs is gener- 
ally described by a double blackbody model, where the emission 
radii (as seen by a distant observer) are typically smaller than 
a few km. The double blackbody, unusual pulse profiles, and 
strongly pulsed fractions have been interpreted in the framework 
of a model, where a pair of thermal, antipodal hot spots of differ- 
ent sizes and temperatures are present at the neutron star surface 
(see Gotthelf et al. 2010, and the case of the CCO RX J0822- 
4300 in the supernova remnant Puppis A). This temperature 
anisotropy can be understood as the result of a large difference 
in intensity between crustal and dipolar components of the mag- 
netic field of the neutron star (Turolla et al. 2011; Shabaltas & 
Lai 2011). As a CCO ages, its host supernova remnant is ex- 



pected to dissipate, whereas the pulsar does not undergo consid- 
erable spin down. If standard cooling is at work (e.g. Page 1998), 
the neutron star can still be observed as a thermal source for the 
subsequent several Myrs. However, intriguing observational ev- 
idence suggests that the CCO CXOU J232327.9+584842 in the 
supernova remnant Cassiopeia A might have experienced en- 
hanced cooling in the past lOyr (Heinke & Ho 2010), which 
has been interpreted as evidence of the onset of neutron su- 
perfluidity in the core of this young neutron star (Page et al. 
2011; Shternin et al. 2011; alternative interpretations exist, e.g. 
Blaschke et al. 201 1). That no "orphan" CCO (i.e. without a su- 
pernova remnant) is recognized among the X-ray emitting INSs 
observed in the solar vicinity, whereas they are expected to rep- 
resent about one-third of the neutron star population (Halpern 
& Gotthelf 2010), might indeed suggest that CCOs cool down 
faster than expected, also perhaps due to accreted light-element 
envelopes once the neutron star is in the photon-cooling stage 
(age > 10 5 yr; Yakovlev & Pethick 2004). 

5.3.2. Relations to peculiar objects and INS populations 

In the following discussion, we consider the thermal X-ray lumi- 
nosities (Lx), emission radii (Rbb), and blackbody temperatures 
(kT) reported in Durant & van Kerkwijk (2006), Mereghetti 
(2008), Ng et al. (201 1, for AXPs), Kaplan (2008, for the M7), 
Halpern & Gotthelf (2010, and references therein; for CCOs), 
Bogdanov et al. (2006, 201 1, for the X-ray properties of MSPs 
in 47Tucanae and M28), Zane et al. (2011, for Calvera), and 
McLaughlin et al. (2007, for RRAT J1819-1458). We note that 
the distance to Calvera is not well-constrained; its spectral prop- 
erties are hence normalised to a distance of 1 kpc. For J 1046, we 
use the corresponding values of Lx = 5x1 32 erg s~ 1 , Rbb = 4 km 
and kT = 135 eV for a neutron star located at d = 2.3 kpc 
(Sect. 2.4). Other pulsar properties are taken from the ATNF 
Pulsar Catalogue (Manchester et al. 2005) 11 . 

1 1 http://www.atnf.csiro.au/people/pulsar/psrcat 
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In Fig. 8, we compare both the X-ray luminosities of these 
neutron stars as a function of emission radius and the emis- 
sion radius as a function of rotation period. The plots show that 
the four main groups - magnetars, CCOs, the M7, and MSPs 
- occupy distinct regions of these diagrams. The Lx - Rbb di- 
agram (Fig. 8; left) shows that the spectral properties of J1046 
are remarkably similar to those of the M7. This resemblance is 
strenghtened by the presence of a broad line in absorption in the 
spectrum of the source (Sect. 5.3.3). On the other hand, the dia- 
gram on the right in Fig. 8 shows a completely different picture: 
the detection of J1046's very rapid rotation is puzzling and at 
odds with its purely thermal, M7-like, spectrum. While all M7 
and the RRAT have spin periods of a few seconds, which are in 
the same range as those observed in magnetars, J 1046 is found 
to be isolated in a region devoid of any other thermally emitting 
INS. 

In general, similar rapidly spinning neutron stars are ex- 
pected to be either young and energetic rotation-powered pulsars 
(with fi dip ~ 10 12 G, E > 10 36 ergs-\ and r ch ~ 10 3 - 10 4 yr) 
or, conversely, old recycled neutron stars with low magnetic 
fields and spin-down luminosities, and old characteristic ages. In 
the absence of any constraint on J1046's spin down, a rotation- 
powered pulsar is a viable, although unlikely, scenario. The ab- 
sence of a supernova remnant or a pulsar wind nebula is indica- 
tive of an age older than a few 10 4 yr and a dipolar magnetic field 
weaker than B dip ~ 7 x 10 11 G. However, the unusually soft and 
thermal energy distribution of the source, its lack of magneto- 
spheric activity, and emission at y-rays are in contrast to other 
rotation-powered pulsars of similar timing properties. 

The likely presence of J 1046 in a star-forming region 
(Sect. 5.2) supports the idea that the neutron star is not old, i.e. it 
was most probably not (fully) recycled in a binary system with 
a low-mass companion that was either evaporated or completely 
accreted (which only happens on timescales of 10 8 yr). If the 
rapid rotation of the source were the result of spin up in a binary 
system, then its companion should have been another massive 
star. In these systems, a shorter period (10 6 - 10 7 yr) of mass 
transfer takes place until the second supernova event disrupts 
the system. J 1046 would then be the first, "mildly-recycled" pul- 
sar ejected from the system, with a spin period and a magnetic 
field in-between those of normal rotation-powered pulsars and 
old MSPs (see e.g. Lorimer et al. 2004; Belczynski et al. 2010). 
Its former companion, which is expected to be a young pulsar, 
could thus still be hidden in the Carina Nebula. We note how- 
ever that the luminosity and spectral properties of J 1046 are still 
atypical relative to those of recycled objects, as illustrated by the 
diagrams in Fig. 8. 

Another intriguing possibility is that of an old CCO, as 
discussed in the case of Calvera (Rutledge et al. 2008; Zane 
et al. 2011). In terms of its spectral properties, Calvera some- 
what differs from J 1046, in that it shows clear evidence of a 
two-temperature thermal model and a non-thermal component 
extending towards higher energies cannot be excluded below 
~ 10% of the source luminosity. As discussed by Zane et al., 
the thermal components possibly originate in two different small 
spots on the stellar surface, a scenario that is consistent with 
Calvera's measurement of a relatively large pulsed fraction. 
These spectral characteristics are those expected from an older 
(i.e. colder and less luminous) and orphan CCO, as also sug- 
gested by the absence of a supernova remnant 12 . Calvera, for a 



As pointed out by Halpern (201 1) and in spite of its spectral proper- 
ties, the upper limit to Calvera's spin down, derived from X-ray timing 
alone, does not exclude a normal rotation-powered pulsar. 



distance of d = 1 kpc, indeed clusters together with the three 
CCOs with known periods in the Rbb - P diagram of Fig. 8. 
Keeping in mind the low statistics and given the still poorly 
known overall characteristics of the population of Galactic anti- 
magnetars, the possibility of an old CCO remains open for 
J 1046, although we note that the overall spectral properties of 
the source do not match those expected for this scenario (e.g. 
double thermal component, size of emission radius, and pulsed 
fraction). 

To more clearly understand the nature of this INS, further X- 
ray observations would be extremely helpful, particularly given 
our ability to constrain the neutron star spin-down rate. For 
this matter, we have been granted new XMM-Newton obser- 
vations for the next period of observation (AOll). With a sec- 
ond epoch, we can constrain spin-down rates of higher than P > 
2.5 x 10 -16 s s- 1 and magnetic fields of B dip > 6 x 10 10 G {2a), 
assuming standard magnetic dipole braking. The corresponding 
lower limit to Z?dip is > 3 Act below the mean field strength of nor- 
mal radio pulsars (Faucher-Giguere & Kaspi 2006) and provides 
sufficient evidence to associate J 1046 with either peculiar, low- 
magnetized objects, such as an old CCO and a mildly-recycled 
pulsar, or the normal rotation-powered population. 

5.3.3. Spectral line in absorption 

New XMM-Newton data on J 1046 may also provide additional 
evidence of the spectral line at e = 1.35 keV (Sect. 2.4). Similar 
spectral features in absorption (in a couple of cases, also har- 
monically spaced) are seen in the spectra of several thermally 
emitting INSs: the CCO IE 1207.4-5209 (Sanwal et al. 2002; 
Mori et al. 2005); the rotating radio transient RRAT J1819-1458 
(McLaughlin et al. 2007); six of the M7 (Haberl 2007); pos- 
sibly 13 Calvera (Zane et al. 2011) and now J1046; these have 
typically been modelled by broad Gaussian absorption lines 
added to a simple blackbody continuum (see, however, the case 
of RX J1308.6+2127; Hambaryan et al. 2011). In the case of 
the M7, the broad features have complex and phase-dependent 
shapes and the line depths often exceed 50%. The presence of 
narrow absorption features, possibly related to either interstel- 
lar or circumstellar highly ionised oxygen, was also reported in 
the co-added RGS spectra of three of the M7 (Hambaryan et al. 
2009; Hohle et al. 2012). 

The spectral absorption features are generally understood 
in terms of the neutron star magnetic field, although its inter- 
pretation is not unique. They can be related to cyclotron tran- 
sitions of either protons (B cyc ~ 10 13 - 10 14 G) or electrons 
(B cyc ~ 10 10 - 10 n G); an alternative explanation would rest 
upon atomic transitions in the outermost layers of the neutron 
star (Lai 2001). 

If the line at energy e - 1.35 keV detected in the spectrum of 
J1046 is interpreted as electron cyclotron absorption similar to 
the case of the CCO IE 1207.4-5209, then the magnetic field of 
the source is expected to be of the order of B cyc — 1.5 x 10 11 G. 
This estimate is based on the assumption that the line is the fun- 
damental and that the gravitational redshift of the neutron star is 
z g = 0.3, according to the relation E cyc = 1.16(B/10 n G)(l + 
Zg) -1 . Alternatively, the feature could be the only first harmonic 
of a fundamental at an energy e ~ 0.65 - 0.7 keV, whose de- 
tection is unfortunately affected by our interpretation of a local 

13 As noted by Zane et al., the contributions of the two thermal com- 
ponents in the spectra of Calvera intersect near 0.65 keV, which is prob- 
ably why the spectral fit can be accommodated by introducing an addi- 
tional feature - either a Gaussian line or an edge - around this energy. 
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oxygen overabundance in the Carina Nebula. In this case, the 
magnetic field would be lower, B = 7 x 10 10 G. These estimates 
can be confronted with the upper limits derived from the pulsar 
spin down. 

5.4. Future prospects 

The knowledge of the pulsar spin down (and its associated quan- 
tities Bfop, r c h, E) is crucial to help us characterise the source 
J1046. Furthermore, a constraint on the period derivative signifi- 
cantly narrows the parameter space of searches for a pulsed y-ray 
signal or even for gravitational waves from the neutron star. 

Gravitational wave observatories have the potential to probe 
several aspects of neutron star physics, in particular providing 
constraints on the state of matter at extreme densities (see e.g. 
Andersson et al. 2011, for a review). Although a detection has 
not yet been reported, the observatories have provided a wealth 
of interesting information on the Galactic neutron star popu- 
lation (e.g. Abbott et al. 2008, 2010; Abadie et al. 2010b,a, 
201 lb, a, 2012; Accadia et al. 2010). The second generation of 
gravitational wave detectors (e.g. Advanced LIGO) might allow 
the first detection (Andersson et al. 201 1). 

Rapidly-spinning and young neutron stars in our Galaxy are 
the best targets to search for a continuous gravitational wave sig- 
nal with ground-based interferometric observatories. The rapid 
rotation of J 1046 implies that it emits continuous gravitational 
waves at a frequency / ~ 110 Hz, where LIGO's S /N is the 
highest (see Abbott et al. 2010, their Fig. 4). In principle, a 
spin-down rate as low as P ~ 8 x 10~ 16 ss~', which may be 
detected with the new XMM-Newton observation planned in 
AOll, places J1046 within the sensitivity of the LIGO's S5 
Scientific Run, considering the whole data set of integration. 
This estimate, of course, simplistically assumes a 100% conver- 
sion of spin-down energy into gravitational waves, regardless of 
how large a neutron star asymmetry is required to power such 
emission. In any case, its rapid rotation, quite young age, and 
relative proximity make this neutron star in the Carina Nebula 
an interesting prospect for gravitational wave detection with 
Advanced LIGO (B. Owen, private communication). 

In spite of many searches and interesting case studies, no 
other thermally emitting INS, presenting exactly the same char- 
acteristics as the M7, has been identified outside the solar 
vicinity to date (see also the case of the high-B radio pulsar 
PSR J0726-2612, which is another likely product of the Gould 
Belt; Speagle et al. 201 1). It is therefore very important for pop- 
ulation studies to understand why there are so many thermally 
emitting sources with similar periods (and presumably ages and 
magnetic fields) in such a small volume. Is this an anomaly 
caused by the Sun's current location close to regions of active 
stellar formation of the Gould Belt or does it really mean that ra- 
dio surveys do miss a large population of INSs, at least as large 
as that of normal radio pulsars? To answer these questions, in- 
vestigations at fainter fluxes as well as population modelling in 
the Galactic scale are needed. 

eROSITA (Extended Rontgen Survey with an Imaging 
Telescope Array 14 ) will extend the ROSAT All-Sky Survey 
(Voges et al. 1999) towards higher energies, with unprecedented 
spectral and angular resolution. The effective area of eROSITA 
is twice that of one XMM-Newton telescope at energies below 
2keV and its sensitivity during the planned all-sky survey will be 
approximately 30 times that of ROSAT. These properties com- 
bined will allow an estimated number of from about 100 to 200 
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new X-ray thermally emitting INSs to be discovered, after four 
years of the eROSITA survey. The above figures are derived us- 
ing the population synthesis model described in Pires (2009) 15 . 
Expectations for the same population based on the model of 
Popov, Boldin, and collaborators predict a more modest num- 
ber, of ~ 30 new discoveries after the planned years of survey. 
Although many of the assumptions considered in our population 
synthesis model differ from those of Popov et al., preventing a 
proper comparison, we estimate that the likely parameters con- 
trolling the discrepancy are the adopted cooling rates and, to a 
lower extent, the model of the interstellar medium distribution. 

6. Conclusions 

Our observational campaign on 2XMM J 104608. 7-594306 has 
confirmed previous expectations and revealed a unique object. 
Most notably, our deep XMM-Newton observation has allowed 
us to unveil the rotation of the neutron star: a very fast spin pe- 
riod of 19 ms is surprisingly at odds with J1046's purely ther- 
mal energy distribution, which is reminiscent of those of the 
Magnificent Seven. Spectral features in absorption have also 
been positively identified. The nature of the source remains open 
to different possible interpretations, a relation to the still poorly 
known class of Galactic anti-magnetars being a favoured sce- 
nario. New XMM-Newton data, granted for the next cycle of 
observations (AOll), will greatly improve the current observa- 
tional picture of this source, and enable us to tightly constrain the 
pulsar spin down. Dedicated radio observations would be simi- 
larly invaluable to help us unveil any truely radio-quiet nature. 
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